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Abstract.  This  paper  discusses  the  short-period  terms  of  nutation  which 
are  included  in  the  recent  theories  of  nutation  for  a  rigid  Earth.  We  argue 
that  these  terms  should  be  expressed  as  terrestrial  perturbations,  that  is 
as  polar  motion,  and  describe  how  such  a  transformation  can  be  accom¬ 
plished  in  practice.  Based  on  the  available  rigid-Earth  amplitudes,  we 
estimate  this  lunisolar  effect  in  polar  motion  for  an  elastic  Earth  with 
liquid  core,  and  compare  it  to  the  oceanic  and  atmospheric  perturbations 
at  similar  frequencies. 


1.  Introduction 

Recent  developments  of  the  rigid-Earth  nutation  theory,  such  as  SMART97 
(Bretagnon  et  ah,  1998),  REN2000  (Souchay  et  ah,  1999),  RDAN97  (Roosbeek, 
1999),  contain  short-period  terms  caused  by  the  interaction  of  the  lunisolar  tidal 
potential  with  those  geopotential  components  which  express  departures  from  ro¬ 
tational  symmetry.  These  are  quasi-diurnal  nutations  corresponding  to  the  6*3,4, 
5*3,4,  T7 4,1 ,  5*4,1  components  of  geopotential,  quasi- semidiurnal  nutations  corre¬ 
sponding  to  the  C 2,2,  S 2,2,  6*3,2,  5*3,2  components,  and  quasi-terdiurnal  nutations 
corresponding  to  the  6*3,3,  5*3,3  components.  Most  of  these  nutations  are  pro¬ 
grade  (counterclockwise),  the  largest  retrograde  amplitudes  being  of  the  order 
of  0.5  pSuS  (microarcsecond).  Therefore,  when  expressed  in  the  terrestrial  frame 
they  become  the  low  frequency  (prograde  and  retrograde),  diurnal  prograde  and 
semidiurnal  prograde  variations,  respectively,  in  polar  motion.  This  part  of  the 
nutation  spectrum  has  received  some  attention  recently  (e.g.  Bretagnon,  1998; 
Bizouard  et  ah,  2000),  nevertheless  it  is  not  included  in  the  nonrigid-Earth  nu¬ 
tation  models. 

The  aim  of  this  paper  is  to  describe  short-period  nutation  and  to  study 
its  importance  for  the  astrometric  determination  of  the  Earth  orientation  pa¬ 
rameters  and  for  their  geophysical  interpretation.  We  start  from  the  kinematical 
considerations  which  describe  how  the  diurnal  and  sub-diurnal  terms  of  nutation 
can  be  expressed  as  polar  motion.  Next,  we  estimate  how  their  amplitudes  are 
modified  by  the  nonrigid-Earth  response.  Finally,  we  compare  this  part  of  the 
lunisolar  effect  on  the  equatorial  component  of  Earth  rotation  with  the  oceanic 
and  atmospheric  perturbations  at  similar  frequencies. 
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2.  Kinematical  relationships 

Current  observational  practice  is  to  express  perturbations  in  Earth  rotation  by 
the  following  set  of  the  Earth  orientation  parameters  { xp ,  yp,  Sip,  Se,  UT1— UTC}, 
where  xp  ,  —  yp  denote  the  terrestrial  coordinates  of  the  Celestial  Ephemeris 
Pole  (CEP),  Sip,  Se  are  the  spatial  offsets  in  longitude  and  obliquity  of  the  CEP 
with  respect  to  its  position  predicted  by  the  precession/nutation  model,  and 
UT1  — UTC  denotes  the  variation  of  Universal  Time  (McCarthy,  1996).  We 
will  consider  polar  motion  as  a  terrestrial  motion  of  the  CEP  and  denote  it 
p  =  xp  —  iyp,  and  nutation  as  a  residual  (that  is  with  respect  to  the  conven¬ 
tional  a  priori  model)  celestial  motion  P  =  X  +  iY  =  sin  s0Sip  +  iSe,  where 
i  =  i/—l  denotes  the  imaginary  unit  and  e0  is  the  mean  obliquity  at  the  refer¬ 
ence  epoch.  For  more  details  concerning  this  parameterization  of  Earth  rotation 
see  Brzezinski  and  Capitaine  (1993). 

2.1.  Nutation  versus  polar  motion 

There  is  a  full  degeneracy  between  polar  motion  and  nutation  parameters:  resid¬ 
ual  nutation  P  (or  its  part)  can  be  equivalently  expressed  as  polar  motion  A p 
using  the  following  first-order  relationship  (Brzezinski  and  Capitaine,  1993) 


A  p=-Pe~te,  (1) 

which  should  be  understood  in  a  sense  that  new  polar  motion  coordinates  are 
p  +  Ap  .  Here  0  =  f lt-\-0o(t)  is  the  rotation  angle  around  the  CEP  axis,  where  12 
is  an  average  angular  speed  of  diurnal  sidereal  rotation  and  60  is  a  quasi- const  ant 
phase  referred  to  t  =  0. 

A  single  quasi-harmonic  term  of  nutation  is  usually  expressed  as 

A  ip  =  ipssm-f(t)  +  ip  c  cos  j(t) ,  (2) 

Ae  =  es  sin  y(t)  +  ec  cos  y(t) ,  (3) 

(e.g.  Souchay  et  al.,  1999),  where  the  nutation  argument 

y(t)  =  kilm  T  k^F  k^D  -\-  k^kl  -\-  k^cp  (4) 

is  integer  combination  of  the  fundamental  arguments:  lm  -  the  mean  anomaly 
of  the  Moon,  ls  -  the  mean  anomaly  of  the  Sun,  F  -  the  difference  between 
the  mean  longitude  of  the  Moon  and  the  mean  longitude  of  its  node,  D  -  the 
difference  between  the  mean  longitudes  of  the  Moon  and  of  the  Sun,  12  -  the 
mean  longitude  of  the  node  of  the  moon,  and  of  the  Greenwich  mean  sidereal 
time  (GMST)  <f>.  This  argument  takes  the  form  y(t)  =  A t  +  7 0(t)  where  A  is  the 
angular  frequency  of  nutation  and  7  0(t)  a  quasi- const  ant  phase.  For  k§  =  0  we 
have  classical  long  periodic  nutation  components,  while  for  /eg  =  1,2,3  diurnal, 
semidiurnal  and  terdiurnal  components  of  nutation,  respectively. 

The  corresponding  perturbation  of  the  CEP  is  an  elliptical  celestial  motion 
which  can  be  expressed  as  a  sum  of  two  quasi-uniform  circular  oscillations  with 
the  same  period  but  in  the  opposite  directions 

P  =  sin  e0  Aip  +  iAs  =  a+  e^+  e1^  +  a_  elv~  e~1^ , 


(5) 
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where  a  >  0  denotes  the  amplitude,  ip  is  the  phase  with  respect  to  the  argument 
7.  If  A  >  0,  which  is  the  case  for  the  short-period  nutation,  the  term  labelled  by 
the  subscript  “+”  describes  the  prograde  (counterclockwise)  motion,  while  that 
with  the  subscript  “  — ”  the  retrograde  (clockwise)  motion.  From  the  comparison 
of  eq.  (5)  to  equations  (2),  (3)  we  derive 

a+  cos  <p+  =  0.5(r/y  sin  e0  +  es),  a+  sin  <p+  =  0.5(— ifis  sin  e0  +  ec),  (6) 
a_  cos  <p_  =  0.5(r/y  sin  e0  —  es),  a_  sin  <p_  =  0.5(V>S  sin  e0  +  ec).  (7) 

For  the  diurnal  and  sub-diurnal  nutation  terms  the  retrograde  amplitudes  are 
much  smaller  than  the  prograde  ones,  the  largest  being  about  0.6  pas,  therefore 
we  will  neglect  them  retaining  in  eq.  (5)  only  the  prograde  term  and  dropping 
the  subscript  “+”  for  simplicity.  After  neglecting  small  differences  between  0 
and  </>,  eq.  (1)  yields  the  following  terrestrial  representation  of  the  nutation 

A p  =  Axp  -iAyp  =  a e^+^'W,  (8) 

with 

7  (f)  =  7(0  ~  41  =  kilm  +  +  k% F  k/±D  +  +  (kg  —  (9) 

Hence,  in  the  conventional  representation  of  nutation,  the  transformation  from 
the  celestial  to  the  terrestrial  representation  is  very  simple:  it  consists  of  de¬ 
creasing  the  coefficient  of  GMST  in  7  by  1  and  reversing  the  phase  ip. 

2.2.  CEP  versus  instantaneous  pole  of  rotation 

Let  u  =  (mi,  m2, 1  +  m3)1  be  the  Earth’s  instantaneous  rotation  vector  ex¬ 
pressed  in  the  terrestrial  reference  system  (TRS),  where  m;  are  small  dimension¬ 
less  quantities  introduced  by  Munk  and  MacDonald  (1960).  Any  perturbation 
of  the  corresponding  pole,  the  instantaneous  pole  of  rotation,  m  =  m\  +  i  m2 
can  be  related  to  the  reported  quantity,  polar  motion  p  of  the  CEP,  using  the 
following  first-order  relationship  derived  by  Brzezinski  (1992) 

P 


or  in  the  frequency  domain  representation 


m( 7 


eicrt 


(T 

O 


Pa  e 


iat 

•> 


(ii) 


where  <7  denotes  the  Earth-referred  angular  frequency  and  ma,  pa  are  ampli¬ 
tudes,  complex  in  general.  See  also  Brzezinski  and  Capitaine  (1993)  for  discus¬ 
sion. 

For  low-frequency  polar  motion  the  differences  between  the  amplitudes  ma 
and  pa  are  small,  but  for  the  prograde  diurnal  and  semidiurnal  variations  the 
ratio  p<Jlmcr  =  (1+  cr/fl)-1  equalling  1/2,  1/3,  respectively. 
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3.  Lunisolar  excitation  of  polar  motion 

Let  us  denote  by  L  =  L\  -\-iL2  the  equatorial  component  of  the  external  gravita¬ 
tional  torque  on  the  Earth,  expressed  in  the  TRS,  which  is  responsible  for  polar 
motion.  Let  us  further  assume  that  it  takes  a  harmonic  form 

L  =  Lael°\  (12) 

where  in  the  present  case  the  frequency  a  is  in  the  vicinity  of  0,  17  or  217. 

We  will  consider  a  very  simple  model  of  the  Earth  consisting  of  the  rotation- 
ally  symmetrical  elastic  mantle  and  the  liquid  core,  with  no  coupling  between 
these  two.  Note,  however,  that  this  model  is  fully  consistent  with  that  applied  to 
estimate  the  corresponding  long-period  and  diurnal  prograde  atmospheric  and 
oceanic  influence  on  polar  motion.  (See  e.g.  Brzezinski  and  Petrov,  2000;  Chao 
et  al. ,  1991;  1996;  Gross  et  al. ,  1998).  The  elastic  yielding  of  the  mantle  is  a 
sum  of  two  different  effects  (Munk  and  MacDonald,  1960): 

•  rotational  deformation  proportional  to  the  polar  motion  amplitude  m  and 
appearing  independently  of  the  physical  nature  of  perturbation,  and 

•  tidal  deformation  associated  with  external  gravitational  forcing  and  propor¬ 
tional  to  the  torque  L. 

After  taking  into  account  these  effects,  the  corresponding  polar  motion  so¬ 
lution  is  obtained  in  the  following  form 

m  =  ma  eirjt ,  (13) 

with  the  (complex)  amplitude  expressed  by 


1  k  d+Q 

1  1  h.  O 


mt r  =  - 


(7  (7e  A mkl  1  -|- 


where  k  denotes  the  tidal  Love  number,  ks  =  0.94  is  the  secular  Love  number, 
em  =  (Cm  —  Am)/Am  stands  for  the  dynamical  ellipticity  of  the  mantle,  Am,  Cm 
are  the  principal  equatorial  and  axial  moments  of  inertia  of  the  unperturbed 
mantle,  and 

1  -  T- 

<?e  =  em  - - J7  (15) 

1  A  jr 

is  the  Chandler  frequency  of  the  mantle.  This  is  eq.  (3-146)  of  Moritz  and 
Mueller  (1987)  in  which  the  principal  moments  A,C  of  the  entire  Earth  have 
been  replaced  by  those  confined  to  the  mantle  alone. 

We  can  compare  this  solution  to  the  corresponding  solution  for  the  rigid 
Earth,  that  is  the  case  k  =  0,  Am  =  A,  Cm  =  C,  by  the  ratio 

~r  =  <l(a)  ?  (16) 


in  which  the  nonrigid  Earth  response  function  q(cr)  is  given  by 


o-or  A 

=  7T 

(7  —  (7  p  J\.  m 


a  1 _ k_  CT  +  fi 

A  r  i  o 


(17) 
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where  ay  is  the  Euler  frequency  of  resonance  for  the  rigid  Earth,  Lra  is  the 

lunisolar  torque  on  the  entire  Earth,  and  La  is  that  on  the  mantle. 

Let  us  first  consider  separately  various  contributions  to  q(cr): 

1.  Elastic  rotational  deformations  proportional  to  m: 

—  Change  of  the  resonant  frequency  from  ay  =  0/304  to  ay  =  0/430.  The 
corresponding  factor  (a  —  ay)/(a  —  ay)  differs  significantly  from  unity  only  at 
low  frequencies. 

—  Factor  [1  +  (k/ks)  em]_1  ~  0.999  which  is  not  important. 

2.  Elastic  tidal  deformations  (proportional  to  L): 

—  Factor  1  —  (k/ks)  (a  +  0)/0  which  equals  about  0.7,  0.4,  0.1  for  a  =  0,  O,  20, 
respectively.  Note  that  for  the  long  periodic  nutations  (a  ~  —O)  this  factor 
is  close  to  1. 

3.  Decoupling  of  the  liquid  core: 

—  Factor  Aj Am  =  1.129. 


Table  1.  Nonrigid  Earth  response  functions:  W81  -  (Wahr,  1981); 
DD1,  DD2  -  (Dehant  and  Defraigne,  1997),  eq.  (4),  eq.  (5),  respec¬ 
tively;  M91  -  (Mathews  et  al. ,  1991).  All  the  functions  have  been 
modified  by  adjusting  the  Chandler  frequency  to  the  observed  value 
of  1  cycle  in  433  days.  Periods  are  in  solar  days  and  amplitudes  in 
microar  cseconds . 


Celestial 

period 

Terrestrial 

period 

Rigid  Nonrigid  Earth  response  coefficients 

ampl.  This  work  W81  DDl  M91 

DD2 

1.07545 

-13.719 

1.23 

0.825 

0.797 

0.803 

0.799 

0.800 

1.03505 

-27.322 

13.60 

0.822 

0.796 

0.801 

0.796 

0.797 

1.00243 

-193.56 

0.63 

0.895 

0.866 

0.873 

0.850 

0.854 

0.99772 

-2190.4 

1.22 

1.070 

1.033 

1.041 

0.981 

0.991 

0.99758 

-3231.5 

7.98 

1.086 

1.048 

1.056 

0.993 

1.004 

0.99727 

OO 

0.74 

1.125 

1.086 

1.094 

1.023 

1.035 

0.99696 

3231.5 

9.61 

1.177 

1.135 

1.144 

1.062 

1.075 

0.99682 

2190.4 

1.62 

1.208 

1.164 

1.174 

1.083 

1.100 

0.99455 

365.24 

0.89 

-1.025 

-0.962 

-0.975 

-0.590 

-0.653 

0.99425 

328.17 

0.53 

-0.264 

-0.237 

-0.243 

-0.019 

-0.056 

0.99216 

193.56 

2.83 

0.516 

0.506 

0.508 

0.566 

0.557 

0.96215 

27.322 

15.27 

0.755 

0.734 

0.738 

0.744 

0.743 

0.92969 

13.719 

1.15 

0.754 

0.734 

0.738 

0.742 

0.742 

0.52743 

1.11951 

2.36 

0.486 

0.495 

0.475 

0.521 

0.522 

i : _ :  „  j. _ _  i  _  j. : _ _ i :  _  _ 

0.49863 

0.99727 

14.31 

0.449 

0.461 

0.436 

0.490 

0.492 

0.34502 

0.52752 

0.06 

0.147 

0.172 

0.093 

0.241 

0.243 

0.33652 

0.50790 

0.05 

0.122 

0.147 

0.063 

0.220 

0.223 

Our  q  is  compared  in  Table  1  to  the  response  functions  used  in  the  recent 
nonrigid-Earth  nutation  theories,  derived  for  more  realistic  multi-layer  structural 
models  of  the  Earth.  These  functions  are  well-constrained  at  the  frequencies 
of  the  long-period  nutations,  which  has  been  verified  by  comparison  with  the 
nutation  amplitudes  observed  by  VLBI,  but  this  is  not  necessarily  the  case  for 
the  short-period  nutations.  In  any  case,  our  q  is  consistent  with  other  estimates 
to  a  few  per  cent.  When  expressed  in  terms  of  the  corresponding  nonrigid-Earth 
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nutation  amplitudes,  the  differences  are  below  0.5  ^as  which  is  a  cut-off  level 
adopted  in  Table  2. 

Let  us  make  the  following  remarks.  First,  when  computing  lunisolar  effects 
on  polar  motion  based  on  the  rigid-Earth  amplitudes  pra  for  the  figure  axis, 
which  is  the  present  case,  the  same  response  coefficient  is  applicable  for  the 
instantaneous  rotation  axis  (Wahr,  1981;  Moritz  and  Mueller,  1987) 

ma  _  pa 

mra  pra  ' 

Secondly,  the  lunisolar  torques  estimates  are  based  on  the  satellite- determi¬ 
ned  geopotential  of  the  entire  Earth,  that  is  to  a  good  accuracy  our  Lra.  When 
adopting  the  assumption  that  the  core  does  not  participate  in  rotational  varia¬ 
tions,  the  coefficient  La/Lra  contains  information  about  the  large-scale  features 
of  the  core-mantle  boundary  shape.  Therefore,  constraining  the  amplitudes  of 
diurnal  and  sub-diurnal  nutations  from  space  geodesy  observations,  which  is  in 
fact  empirically  determining  the  ratio  La/Lra ,  would  provide  important  infor¬ 
mation  about  the  internal  constitution  of  the  Earth. 


4.  Results  and  comparison  with  oceanic  and  atmospheric  effects 


Now,  we  will  consider  in  detail  diurnal  and  semidiurnal  nutation  components 
expressed  as  low  frequency  and  prograde  diurnal  variations  in  polar  motion, 
respectively.  We  will  disregard  terdiurnal  nutation  components  because  their 
amplitudes  are  below  0.1  pas  for  a  rigid  Earth,  and  according  to  Table  1  this 
number  is  further  decreased  about  8  times  for  an  elastic  Earth  with  liquid  core. 


4.1.  Prograde  diurnal  polar  motion 

Prograde  diurnal  polar  motion  corresponding  to  semidiurnal  nutation,  is  caused 
by  the  interaction  of  the  lunisolar  tidal  potential  with  the  Cn<2,  SU) 2  (n  >  2)  terms 
of  geopotential,  where  only  the  n  =  2  contributions  are  larger  than  0.15  pas. 
This  phenomenon  was  described  and  estimated  for  the  first  time  by  Chao  et  al. 
(1991)  as  the  so-called  “prograde  diurnal  libration”  in  polar  motion.  It  is  driven 
by  the  lunisolar  torque  taking  the  following  form  (see  ibid.) 


L 


3  iGM(B-A) 
4  r3 


sin  20  e 


i(A-A0) 

1 


(18) 


where  G  denotes  gravitational  constant,  M  is  the  disturbing  point  mass,  r,  0,  A 
are  its  geocentric  radius,  co-latitude  and  longitude,  respectively,  B  >  A  are  the 
principal  equatorial  moments  of  the  Earth’s  inertia,  and  A0  denotes  the  longitude 
of  the  A-axis. 

Our  equations  (14)  and  (11)  with  a  =  fl  yield  immediately  the  following 
solution 

p=  -0.168  sin 20  e-*(A-Ao),  (19) 

A  r3Q2  v  7 

which  describes  prograde  a  nearly  diurnal  oscillation  because  A  ~  —  flf+const. 
Note,  however,  that  Chao  et  al.  (1991;  1996)  estimated  the  coefficient  0.36 
instead  of  our  0.168,  which  is  also  reflected  in  the  individual  amplitudes  given  in 
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Table  2.  This  could  be  either  due  to  the  fact  that  their  solution  is  for  m  instead 
of  the  observed  p  (see  Sec.  2.2),  or  that  they  neglected  the  influence  of  the  elastic 
tidal  deformation  on  q  (see  Sec.  3),  both  effects  giving  a  factor  of  about  2. 

This  prograde  diurnal  variation  in  polar  motion  is  proportional  to  the  equa¬ 
torial  dynamical  flattening  ( B  —  A)/A  of  the  mantle,  which  can  differ  from  that 
of  the  entire  Earth  by  as  much  as  40%  (Herring  and  Dong,  1994).  If  it  was  pos¬ 
sible  to  estimate  its  amplitude  from  the  observation  of  polar  motion,  we  could 
constrain  in  this  way  the  62,27  <$2,2  components  of  the  core-mantle  boundary. 


Table  2.  Largest  long  periodic  and  diurnal  prograde  variations  in  po¬ 
lar  motion  (amplitude>  0.5  pas)  corresponding  to  diurnal  and  semidi¬ 
urnal  nutations  from  the  REN2000  tables  (Souchay  et  al. ,  1999).  The 
last  column  shows  estimates  of  Chao  et  al.  (1991,  1996)  (only  of  the 
amplitude  because  the  phase  is  the  same  as  in  our  estimate).  Periods 
are  in  solar  days  and  amplitudes  in  microarcseconds. 


Fundamental  arguments 

Tidal 

Terrestr. 

Amplitudes 

Chao 

l  m 

h 

F 

D 

Q 

cf)  code 

period 

rigid 

elastic 

Phase 

et  al. 

-1 

0 

-1 

0 

-1 

0 

-13.719 

1.23 

1.02 

-83° 

0 

0 

-1 

0 

0 

0 

-27.212 

2.16 

1.77 

-83° 

0 

0 

-1 

0 

-1 

0 

-27.322 

13.60 

11.18 

-83° 

0 

0 

-1 

0 

-2 

0 

-27.432 

0.70 

0.57 

97° 

-1 

0 

-1 

2 

-1 

0 

-193.56 

0.63 

0.56 

-83° 

1 

0 

-1 

0 

0 

0 

-2190.4 

1.22 

1.31 

-83° 

1 

0 

-1 

0 

-1 

0 

-3231.5 

7.98 

8.66 

-83° 

0 

0 

0 

0 

0 

0 

00 

0.74 

0.83 

-140° 

-1 

0 

1 

0 

1 

0 

3231.5 

9.61 

11.31 

97° 

-1 

0 

1 

0 

0 

0 

2190.4 

1.62 

1.95 

97° 

0 

0 

1 

-1 

1 

0 

365.24 

0.89 

0.91 

-83° 

-1 

1 

1 

0 

1 

0 

328.17 

0.53 

0.14 

-83° 

1 

0 

1 

-2 

1 

0 

193.56 

2.83 

1.46 

-83° 

0 

0 

1 

0 

2 

0 

27.432 

0.99 

0.75 

97° 

0 

0 

1 

0 

1 

0 

27.322 

15.27 

11.52 

-83° 

0 

0 

1 

0 

0 

0 

27.212 

2.40 

1.81 

-83° 

1 

0 

1 

0 

1 

0 

13.719 

1.15 

0.87 

-83° 

Sum  of  all  amplitudes 

70.33 

62.76 

-1 

0 

-2 

0 

-2 

1  Q 1 

1.11951 

2.36 

1.14 

-120° 

0.8 

0 

0 

-2 

0 

-1 

1 

1.07598 

2.17 

1.03 

-120° 

0 

0 

-2 

0 

-2 

1  Oi 

1.07581 

11.56 

5.47 

-120° 

11.3 

-1 

0 

0 

0 

0 

1  Ml 

1.03472 

0.83 

0.38 

60° 

0.8 

0 

0 

-2 

2 

-2 

1  Pi 

1.00275 

4.77 

2.15 

-120° 

4.7 

0 

0 

0 

0 

0 

1  K 1 

0.99727 

14.31 

6.43 

60° 

14.6 

0 

0 

0 

0 

-1 

1 

0.99712 

1.94 

0.87 

60° 

1 

0 

0 

0 

1 

1 

0.96257 

0.78 

0.34 

60° 

Sum  of  all  amplitudes 

43.45 

19.99 

Eq.  (19)  can  be  further  refined  by  expressing  r,  0,  A  as  functions  of  time. 
Instead,  we  used  the  REN2000  nutation  tables  for  a  rigid  Earth  (Souchay  et  al., 
1999)  together  with  the  nonrigid-Earth  response  function  given  by  eq.  (17)  and 
assuming  La  =  Lra.  Note  that  this  choice  of  the  rigid-Earth  nutation  theory  has 
no  significant  influence  upon  our  estimate  because  according  to  Bizouard  et  al. 
(2000)  the  agreement  between  the  3  available  theories  mentioned  in  Sec.  1  is  at 
the  level  of  1  pas  in  the  time  domain.  From  the  estimated  amplitudes  of  the 
prograde  diurnal  polar  motion,  given  in  Table  2,  we  can  conclude  that: 
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•  The  maximum  peak-to-peak  variation  for  the  nonrigid  Earth  is  about  40  ^as 
with  the  main  contributions  at  frequencies  of  the  principal  diurnal  tides  K\,  0\ 
and  Pi. 

•  Comparison  with  the  ocean-tide  contributions  (Table  3)  shows  that  prograde 
diurnal  librations  are  systematically  about  25  times  smaller  and  differ  in  phase 
by  about  +90°. 

•  Comparison  with  the  prograde  diurnal  atmospheric  tides  in  polar  motion  (Ta¬ 
ble  3)  shows  similar  order  of  magnitude  but  different  power  distribution  among 
diurnal  frequencies. 

•  Our  estimate  confirms,  and  even  strengthens,  the  conclusion  of  Chao  et  al. 
(1996)  that  present  measurement  and  modeling  accuracies  are  inadequate  for 
the  diurnal  libration  in  polar  motion  to  close  the  budget,  even  for  the  largest 
tidal  waves. 


Table  3.  Prograde  diurnal  and  semidiurnal  polar  motion  excited  by 
the  oceanic  tides  (Chao  et  al.,  1996)  and  by  the  atmospheric  tides 
(Brzezinski  and  Petrov,  2000).  Periods  are  in  solar  days  and  amplitudes 
in  microarcseconds.  Note  that  the  phase  is  referred  to  GMST,  which 
differs  from  the  usual  convention  in  representing  tidal  contributions. 


Fundamental  arguments 

Tidal 

Terrestr. 

Oceanic 

Atmospheric 

l  m 

Is 

F 

D 

Q 

4> 

code 

period 

ampl. 

phase 

ampl. 

phase 

-l 

0 

-2 

0 

-2 

l 

Qi 

1.1195 

27 

167° 

0 

0 

-2 

0 

-2 

l 

O  i 

1.0758 

142 

160° 

0 

0 

-2 

2 

-2 

l 

Pi 

1.0027 

56 

153° 

1.3 

-133° 

0 

-1 

0 

0 

0 

l 

Si 

1.0000 

7.1 

47° 

0 

0 

0 

0 

0 

l 

K  i 

0.9973 

171 

-27° 

1.6 

117° 

0 

1 

0 

0 

0 

l 

tpi 

0.9946 

0.7 

-135° 

-1 

0 

-2 

0 

-2 

2 

n2 

0.5274 

17 

-45° 

0 

0 

-2 

0 

-2 

2 

m2 

0.5175 

75 

-64° 

0 

0 

-2 

2 

-2 

2 

s2 

0.5000 

29 

-94° 

2.9 

-102° 

0 

0 

0 

0 

0 

2 

to 

0.4986 

7 

-84° 

4.2.  Low  frequency  lunisolar  polar  motion 

Low  frequency  lunisolar  polar  motion  corresponding  to  diurnal  nutation,  is 
caused  by  the  interaction  of  the  lunisolar  tidal  potential  with  the  Cn<i,  Sn< \ 
(n  >  3)  terms  of  geopotential,  where  only  the  n  =  3  contributions  are  larger 
than  0.75  //as.  Detailed  expansion  (Table  2)  shows  that  the  maximum  peak-to- 
peak  variation  for  the  nonrigid  Earth  can  reach  120  //as,  which  is  about  3  times 
larger  than  prograde  diurnal  libration  and  significant  in  view  of  the  current  ac¬ 
curacy  of  monitoring  polar  motion.  Individual  contributions  are  not  exactly  at 
the  frequencies  of  important  low  frequency  tides  and  are  well  separated  from  the 
much  larger  atmospheric  contributions  at  seasonal  frequencies.  The  spectrum  is 
approximately  symmetrical  with  respect  to  zero  frequency.  Therefore  the  waves 
are  highly  elliptical.  It  can  be  seen  when  writing  the  largest  waves  explicitly: 

27.32  d:  e_*83°[22.7/nas  cos(F  +  Cl)  +  0.3/uas  sin(_F-|-D)], 

27.21  d:  e_*83°3.6/nas  cos  F, 

8.85  y:  e_*83°[20.0/nas  cos(_F-|-D  —  lm)  +  2.6/tras  sin(_F-|-0  —  lm)\. 
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Comparison  with  the  ocean-tide  contribution  by  the  Mm  wave  done  in  Table  4 
shows  that  the  direct  lunisolar  effect  at  a  monthly  period  is  only  about  4  to  5 
times  smaller  and  is  of  a  size  similar  to  the  current  accuracy  of  determinations 
of  the  amplitudes  from  the  geodetic  observations. 


Table  4.  Comparison  of  the  direct  lunisolar  effect  in  polar  motion  at  a 
monthly  period  with  the  ocean-tide  contribution  and  with  the  observed 
polar  motion  with  atmospheric  contribution  (PM/AAM)  subtracted, 
reported  by  Gross  et  al.  (1998).  Amplitudes  are  in  microarcseconds. 


Retrograde 
ampl.  phase 

Prograde 
ampl.  phase 

This  work  (27.43d) 

0.6 

97° 

0.8 

97° 

This  work  (27.32d) 

rr.2 

-83° 

77.5 

-83° 

This  work  (27.21d) 

r.s 

-83° 

7.8 

-83° 

Ocean  Mm  tide  (27.55d) 
Brosche  /  Seiler  /  Gross 
Dickman’s  model 

55 

28° 

53 

-706° 

r7 

-7° 

32 

-44° 

PM/AAM  data  Mm 

59 

-83° 

44 

-85° 

irs 

±77° 

±20 

±26° 

5.  Conclusions 

Diurnal  and  sub-diurnal  terms  of  nutation  that  are  caused  by  the  interaction  of 
the  lunisolar  tidal  potential  with  those  geopotential  components  which  express 
departures  from  the  rotational  symmetry,  are  of  geophysical  interest  because 
their  amplitudes  depend  on  the  large-scale  features  of  the  core-mantle  boundary. 
We  developed  an  analytical  model  of  the  short-period  nutation  for  a  nonrigid 
Earth,  based  on  the  rigid-Earth  nutation  theory  REN2000  assuming  that  the 
Earth  consists  of  an  elastic  mantle  and  a  liquid  core  with  no  coupling  between 
these  media.  This  model  is  expressed  in  the  terrestrial  frame,  that  is  as  polar 
motion,  which  is  consistent  with  the  treatment  of  the  oceanic  and  atmospheric 
effects  appearing  at  similar  frequencies.  When  taking  into  account  the  size  of 
various  components,  we  can  conclude  the  following: 

•  The  prograde  semidiurnal  polar  motion  corresponding  to  the  terdiurnal  nuta¬ 
tion  is  of  the  order  of  0.01  fj, as,  and  therefore  completely  negligible. 

•  The  prograde  diurnal  polar  motion,  sometimes  referred  to  as  diurnal  libration, 
corresponding  to  the  semidiurnal  nutation  has  a  maximum  peak-to-peak  size  of 
about  40  /was.  There  is  a  little  chance  that  in  the  near  future  this  motion  can  be 
observationally  separated  from  the  much  larger  ocean-tide  influences  appearing 
at  the  same  frequencies. 

•  The  long-period  polar  motion  corresponding  to  the  diurnal  nutation,  with 
main  constituents  at  27.3  days  and  8.85  years,  is  of  the  maximum  peak-to-peak 
size  of  about  120  /ras.  It  is  significant  in  view  of  the  current  accuracy  of  space 
geodetic  observations,  and  can  be  separated  from  the  atmospheric  and  oceanic 
perturbations.  Therefore  it  should  be  accounted  for  in  the  excitation  balance  of 
polar  motion. 
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